BBGKY Hierarchy Underlying Many Particle Quantum Mechanics 
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Recently, the one particle quantum mechanics has been obtained in the framework of an en- 
tirely classical subquantum kinetics. In the present paper we argue that, within the same scheme 
and without any additional assumption, it is possible to obtain also the n-particle non relativistic 
quantum mechanics. The main goal of the present effort is to show that the classical BBGKY hier- 
archical equation, for the n-particle reduced distribution function, is the ancestor of the n-particle 
Schrodinger equation. On the other hand we show that within the scenario of the subquantum 
structure of quantum particle, the Fisher information measure emerges naturally in quantum me- 
chanics. 
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In ref . , it is shown that quantum mechanics can be 
obtained in a self consistent scheme of an entirely classi- 
cal many body physics. Within this framework the mean 
features of quantum mechanics (i.e. the probabilistic na- 
ture of the quantum description, the quantum potential, 
the Schrodinger equation, the quantum operators, the 
Heisenberg uncertainty principle) were obtained starting 
from a classical subquantum kinetics in the phase space 
and without invoking any additional principle. Further- 
more the fundamental constant h emerges naturally as 
an integration constant and represents a free parameter 
for the theory. Within this theory, the astonishing sce- 
nario of a subquantum structure for the quantum parti- 
cle emerges spontaneously. The quantum particle turns 
out to have an internal structure and a spatial dispersion 
and appears to be composed by N identical point like 
and interacting subquantum objects, the monads, obey- 
ing to the laws of Newtonian physics. The theory can be 
viewed as describing a mechanism which permits to con- 
struct the quantum particle starting from its constituents 
and then quantum mechanics appears to describe a phys- 
ical reality. 

The statistical ensemble of these N monads is de- 
scribed in the phase space through the distribution func- 
tion / = f(t,x,v) whose time evolution is governed by 
the kinetic equation Cf — C(f), being C the Liouville 
operator and C(f) the collision integral that takes into 
account the intermonad interactions. We don't make any 
assumption about the structure of the collision integral 
and the nature of the interaction between the monads 
except that, during the point-like collisions the monad 
number, momentum and energy are conserved. The pro- 
jection of the above phase space kinetics into the physical 
space produces a hydrodynamics which naturally leads 
to the one particle spinless quantum mechanics when the 
external force field is conservative. 

Recently in ref. Q , it is shown that the above kinetic 



approach can be easily adopted to treat the quantum 
particle with spin in an external electromagnetic field. 
Of course in this way the spinning quantum mechanics 
can be obtained. 

The energy conservation during the collisions consid- 
ered in the present scheme plays a very important role 
and could be also taken into account in the stochastic 
models of quantum mechanics as it has been recently ob- 
served in ref. 0. 

A question which naturally arises at this point is if 
it is possible to extend the kinetic formalism of ref. pj 
in order to treat non relativistic many particle quantum 
systems. This point is neglected in the literature where 
almost the totality of works devoted to the derivation and 
interpretation of quantum mechanics concern the single 
particle. 

The main goal of the present paper is to show that 
it is possible to obtain, within the theory developed in 
ref. [lj and without any additional assumption, also the 
n-particle quantum mechanics starting from the classi- 
cal kinetics governed in the phase space by the BBGKY 
hierarchy. 

BBGKY Hierarchy: Let us consider a statistical en- 
semble constituted by J\f identical and indistinguished 
point-like particles (monads) obeying the laws of classi- 
cal physics. It is well known that the n-particle reduced 
distribution function /„ = f n (t, x\, V\, X2, v%, x n , v n ) 
obeys the following BBGKY hierarchical evolution equa- 
tion m 



fj, dv. t 



i 



being Fi = —dVi/dxi and 

n 

Vi = V^(x i ) + -J2^ n) i\^-^\) ■ (2) 
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Ci(/n+i) are the collision integrals which can be ex- 
pressed in term of the function f n +i- The assumption 
that during the point-like collisions the monad number, 
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momentum and energy are conserved, implies that the The density of energy is given by 
three functions gi(vi) = 1, 5 2 («i) = «i and g 3 (vi) = vj 
are the collision invariants of Cj(/ n +i) 



1 2 
E = - fx <v > v 



V=\pu 2 



V 



gk(vi) Cj(/ n +i) d 3 v; = ; fc = 1, 2, 3 . (3) being e > the density of the internal energy. 



Let us consider the 6n-dimensional phase 
space (x,v) being x = (xi, X2, x n ) and 
v = (v%,V2, —,v n ). We introduce the two gradient 
operators d/dx = (d/dxi, d/dx 2 , d/dx n ) and 
d/dv = (d/dvi, d/dv2, d/dv n ). Then we define the 
3n-dimensional force F = (F\, F2, F n ) which can be 
derived from a potential 



dV 
dx 



; V = £> 



(4) 



The BBGKY hierarchical equation (JJ, governing the 
evolution of the n-particle reduced distribution function 
fn = fn(t, x, v), can be written in the following compact 
form 

d ^ ,F df n , . 

^ + ^^ + 7T'^r = c(/ " +l) ■ (5) 

The collision integral 

71 

G(/«+i) = x;a(/n+i) , (6) 

i=l 

admits the three collision invariants gi(v) = 1, (72(f) = w 
and 53(f) = f 2 satisfying the conditions 



g k (v)C(f n+1 )d an v = ; fc = l,2,3 



(7) 



as one can verify immediately starting from Eq. Q . 

Hydrodynamics: We consider now the projection of the 
kinetics from the 6n-dimensional phase space into the 3n- 
dimensional space of the coordinates x, where the system 
is described through the distribution function p n (t, x) = 
J f n (t,x,v) d 3n v. After recalling that the total value of 
a given physical quantity with density A = A(t,x,v) 
can be calculated as J A f n d Sn vd 3n x we define the mean 
value of A in the point x through 



< A(t,x,v) > v 



J A(t,x,v)f n {t,x,v)d 3n v 
J f n (t,x,v)d 3n v 



(8) 



This mean value represents the density of A in the space 
of the coordinates and permits to write the total value 
of A as J p n < A > v d 3n x. We define the densities in 
the space of the coordinates for some quantities used in 
following. The densities of current, of stress tensor and 
of heat flux vector are defined, respectively, by 



u = < v > v , 

<Ty = /i < (Vi - Ui)(Vj - Uj) > v , 

hi = -fi <\v — u\ 2 (vi - Ui) > v . 



(9) 
(10) 

(11) 
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■ Jti (< V >v 



< v >t 



(12) 



(13) 



Using the standard procedure consisting in multiply- 
ing Eq. J5J by the three collision invariants gk(v) and 
integrating with respect to v, the three following hydro- 
dynamic equations can be obtained 0, Q : 



-dT + dx-- (p " u 

Du 

A* 







Dt 



= F + F 



(14) 
(15) 
(16) 



In the above equations D/ Dt — d/dt + u ■ d/dx is the 
total time or substantial or Lagrangian derivative, JF is 
the stress force 



}_d_ 

Pn 9Xj 



(Pn <7ij) 



(17) 



and Si = Eui+UijUj + hi is the energy flux density vector. 

Eq.s 114|) - (|16|) define the most general hydrodynamics 
of the particle system which behaves in the physical space 
as a fluid. In particular Eq. (|16|l imposes the conserva- 
tion of the total energy H of the system 



H 



Ep n d 3n x 



(18) 



After recalling that the external force F is conserva- 
tive, from Eq. (|15|l it follows immediately that the ir- 
rotationality of u implies the conservativity of J- and 
viceversa, namely 



u — 



IdS 
/i dx 
dW 
dx 



(19) 
(20) 



The condition (|19|1 imposes that the system here consid- 
ered is spinless and reduces the vector equation l|15|) into 
a scalar one 



dS 



1 



dt 2p, \dx 



dS 



+W+ V = 



(21) 



Quantum potential: From the two expressions of the 
force J- given by Eq.s l|17|) and l|20|l it follows immediately 



dW 

dx 3 



da 



dxi 



dxi. 



(22) 



with £ = In p n . Note that Eq. (12211 descends from the 
requirement that the force J- is conservative and can be 



3 



viewed as a condition constraining the forms of W and 
<Tjfc which result to be two Junctionals of the field £. The 
solutions of the latter equation are couples of W and ajk 
and in principle there can exist more than one solution. 
Even though it can appear that Eq. (|22|1 contains consid- 
erable degrees of freedom, in the following we will show 
that the particular structure of this equation restricts 
strongly the number of its solutions. 

We suppose now that just as the left hand side, both 
the first and the second term in the right hand side in 
Eq. IL'l'li are curl free and will take the form d(...)/dxj. 
This requirement for the first term can be satisfied in 
two different ways. The first and most natural way is 
related to the choice ajk = aSjk- Alternatively we can 
pose ajk — dak/dxj and after recalling the symmetry 
property ajk — akj we have that ak = da/dxk- Then 
the most general form of the density of stress tensor is 



ajk = aSjk + 



d 2 a 
dxjdxk 



(23) 



being a and a two scalar functionals depending on £. 
After taking into account the identity 



d 2 a <9£ 
dxjdxk dxk 



d ( da <9£ \ da <9 2 £ 



(24) 



dxj \dxk dxkj dxk dxkdxj 

Ea. (l22|l can be written as 

d / , d 2 a da d£ 
[W - a — 

dxj \ dxkdxk dxk dxk 

da d 2 £ 
dxj dxk dxkdxj 

Clearly the two terms in the right hand side of Eq. 125|) 
must be curl free. These requirements impose the two 
following conditions 



(25) 



da d£, 

dxi dx m 
d 2 £ 

dxidxk dxkdx m 



d 2 a 



da d£ 

dx m dxi 
d 2 a 



d 2 i 



dxmdxk dxkdxi 



(26) 
(27) 



It is trivial to verify that from Eq. (|26[1 it results that 
a = er(£) is an arbitrary function of £ while from Eq. 
(|27|l one can obtain that a — c£ with c a real arbitrary 
constant. With these positions Eq. I|25|l becomes 



d 

dxj 



W 



adi 



d£ dZ 



dxkdxk 2 dxk dxi 



d 2 i 







. (28) 

From this last equation we can obtain the expression of 
W while Eq. (|23|) gives the expression of <7y. Finally, 
we can write the most general solution of Eq. (1221) under 
the form 



a jk 



d 2 lnp» 

dxj dxk 

d 2 
dx 2 



5jk rJ f 



ln P„ 



dx 



dp n , (29) 
U( Pn ) . (30) 



The above expressions of ajk and W define a family of 
solutions of Eq. l|2^1) and represent the hydrodynamic 
constitutive equations for the system. Remark that these 
equations have been enforced exclusively from the fact 
that the internal stress forces are conservative. 

In the above constitutive equations the terms contain- 
ing U describe the standard classical Eulerian fluid. For 
instance the nonlinear potential U{p n ) — ap n appearing 
in the so called cubic non-linear Schrodinger equation is 
originated from the stress tensor ajk = \ ap n 8jk- Anal- 
ogously the logarithmic non-linearity U{p n ) = bhip n in 
the Schrodinger equation is originated from a constant 
stress tensor ajk = b Sjk ■ 

In the following, we pose U = and focalize our at- 
tention to the term of Eqs Ij29(l and (|3t)|> . proportional 
to the arbitrary constant c, which is absent in classi- 
cal kinetics. From the definition of the total internal 
energy, it results that H = J d 3n xp n e > 0. Starting 
from Eq. (|29|l and performing an integration by parts, 
we obtain H = —cf d 3n x {dp n /dx) 2 / p n . Then we pose 
c = —rf/Ap, in order to put H > 0. The real positive 
constant rj remains a free parameter of the theory. 

We introduce the normalized n-particle reduced dis- 
tribution function g n = p n j N being N — J p n d 3n x, so 
that J g n d 3n x = I. Then we define the mass m = Np 
the internal potential W = NW and set Ti = Nrj. From 
(|30|l follows the expression of W 



W 



,-1/2 



2m 



9 2 1/2 

dx 2 



(31) 



We recall that the potential W given by Eq. (1.'il[) is the 
quantum potential (Madelug 1926) and it is remark- 
able that here it has been obtained in the framework of 
an entirely classical subquantum kinetics. Note that the 
subquantum monad structure of the system generates the 
stress forces and then the potential W. Clearly when this 
structure is suppressed the quantum potential vanishes. 
Furthermore the presence of the collision integrals C'i in 
the BBGKY kinetic Eq. which take into account the 
monad interactions are necessary for the consistency of 
the theory. 

Many particle Schrodinger equation: We define the ex- 
ternal potential V — NV, and set S = NS. Now the 
quantum fluid is described completely by the two scalar 
fields g n and S whose evolution equations are l|32|) and 
(O, respectively. 



dt 



d_ 

dx 

dS_ _L_ fdS_ 
dt 2m V dx 







" m dx y 

2 

+W+V =0 



(32) 
(33) 



Alternatively, we can describe this fluid by means of the 
complex field <£„ = q\[ 2 exp(iS/h) with J |^„| 2 d 3n a; = I, 
whose evolution equation 



ih 



djn 

dt 



h 2 d 2 * n 

2m dx 2 



(34) 
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can be obtained directly by combining Eq.s i|32|) and (|33|l . 
This later equation is the n-particle Schrodinger equation 
which, after recalling V = J2 i ^ with Vi = NVi, assumes 
the form 



from which it results / = ^ . Ii and then TL = J2i % ■ An 
analogous property holds also for the quantum potential 
W = Wi being 



1T1 



at 



E 

i=l 



2m dx 2 



(35) 



The total energy H given by Eq. (|18|l results to be the 
sum of two terms: H — + TL. The first term 



H (d) 



i=l 



+ > /d^are^x), (36) 



corresponds to the energy of n interacting classical parti- 
cles while the second term TL = J ep n d 3n x is originated 
by the internal structure of the n bodies and is given by 



TL 



8m 



j3n. 



1 

Qn 



dx 



8m 



(37) 



Eq. (|37|) hints on the existence of a tight link be- 
tween quantum mechanics and Fisher information the- 
ory. It is remarkable that TL results to be proportional to 
the Fisher information measure / (Fisher 1922) which 
describes an important feature of quantum mechanics 
[U EH EH E3- Once again, within the scenario of the 
subquantum structure of quantum particle, another ba- 
sic physical concept, namely /, emerges naturally. In the 
present context / has a very transparent meaning be- 
ing simply a measure of the internal energy of quantum 
system. The fact that TL is originated from the internal 
structure of the n-bodies becomes more transparent if we 
observe that TL can be written also as 



-2 

8m 



(38) 



Wi = -—Q. 

2m 



-1/2 



02 1/2 

dx 2 , 



(39) 



Of course starting from H = +TL and taking into ac- 
count Eq.s (|35|l and (|36|l one obtains immediately the well 
known expression of the Hamiltonian of the rt-particle 
quantum mechanics 



H 



-73n „ 



2m 



dx 



V\^ n \A=Y,H l (40) 



i=i 



In the present picture the mass of the n-body quantum 
system is M — nm. After recalling that m = Nfi follows 
immediately that the system is composed by N = nN 
monads. 

Conclusions: The approach used here to obtain the 
Schrodinger equation for a system of n quantum par- 
ticles starting from the BBGKY hierarchy describing a 
statistical system of M = nN subquantum monads, sug- 
gests that quantum mechanics can be viewed as a hidden 
variables theory. This approach can be used easily also to 
construct the non linear quantum mechanics [it is enough 
to consider U ^ in Eqs i(2'9"|) . lj30"|l ] and then for study the 
non linear effects originated for instance from the particle 
bosonic or fermionic nature 0, 0] . 
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